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An underwater acoustic sensor array spanning ^^1500 km'' is used to search for cosmic-ray neu- 
trinos of ultra-high energies (UHE, E^ylO^^ eV). Approximately 328 million triggers accumulated 
over an integrated 130 days of data taking are analysed. The sensitivity of the experiment is de- 
termined from a Monte Carlo simulation of the array using recorded noise conditions and expected 
waveforms. Two events are found to have properties compatible with showers in the energy range 
10^* eV<Esh<5xlO^'^ eV and 10^^ eV<Esh<5xW^'^ eV. Since the understanding of impulsive back- 
grounds is limited, a flux upper limit is set providing the most sensitive limit on UHE neutrinos 
using the acoustic technique. 

PACS numbers: 95.85.Ry, 43.60.Bf, 14.60.Lm 
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I. INTRODUCTION 



In the last few decades considerable attention has been 
devoted to understanding the origin, acceleration mech- 
anisms, flux, and composition of the highest energy cos- 
mic rays. Neutrinos are thought to play an important 
role at ultra-high energies (UHE, E>10^® eV). Different 
mechanisms are being explored to probe the extremely 
small neutrino flux (<^ 1 km~^ yr~^) expected in this 
regime. In the last decade, feasibility studies for acoustic 
UHE neutrino detectors have been initiated in large nat- 
ural bodies of water, ice and salt. The advantage to this 
technique is the possibility of building very large arrays 
(thousands of km'^) with sparse microphones, thanks to 
the large attenuation length of sound at the appropri- 
ate frequencies in these media. The Study of Acoustic 
Ultra-high Energy Neutrino Detection (SAUND) phase 
II is the first experiment to read out hydrophones un- 
dersea for the purpose of detecting UHE neutrinos using 
such an expansive array (^^1500 km'^), and follows the 
first phase [l| where ^15 km'^ were read out at the same 
site. 



Acoustic radiation is produced from the volume ex- 
pansion caused by the heating of the medium where the 
neutrino interacts and produces a shower. Theoretical 
models ^ and experimental measurements Q show that 
the signature of this event is a single bipolar pulse, with 
acoustic energy concentrated around 10 kHz. Noise con- 
ditions measured by SAUND II have been studied in de- 
tail [J] . The analysis revealed that although the ambient 
noise correlates well with surface wind speeds, at higher 
frequencies the spectrum rolled off much faster than the 
expected power law due to the large depth of the hy- 
drophones. The expected noise profile was parametrized 
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and a good agreement was observed when compared to 
measured profiles. 

Full data processing of pulse shapes, analysis of multi- 
hydrophone coincidences, and other techniques intended 
to reject background are reported here, along with a UHE 
neutrino flux limit set by a detailed simulation of the full 
experimental procedure. 



II. THE SAUND II EXPERIMENT 

The second phase of SAUND (SAUND II) employs a 
large hydrophone array at the US Navy's Atlantic Un- 
dersea Test and Evaluation Center (AUTEC) located 
at the Tongue of the Ocean (TOTO) in the Bahamas. 
SAUND II is currently the largest test for the feasibility 
of acoustic ultra-high-energy neutrino detection. This 
program follows a general study of the expected perfor- 
mance 5], and a first experimental phase (SAUND I) 
using seven hydrophones [l|. Since then, the US Navy 
has upgraded the hydrophones and the readout system 
of the array. SAUND II uses 49 of these hydrophones 
with digitized signals transmitted to shore over optical 
fibers. The array spans an area of ^ 20 km x 50 km 
with spacing of 3 to 5 km. Hydrophones are mounted 
5.2 m above the ocean floor, at depths between 1340 and 
1880 m. Fig. [T] shows the configuration and the topogra- 
phy of the ocean fioor. Hydrophones are omnidirectional 
with a flat response (within 5 dB) at the frequencies con- 
sidered. The gains of the 49 channels coincide to within 
1 dB. Analog signals are regenerated at the shore station 
from the digital data (by the US Navy for their backwards 
compatibility) and fed to the SAUND II data acquisition 
(DAQ) system that re-digitizes them at 156 kHz. Since 
low frequencies are not relevant for SAUND II, a hard- 
ware high-pass RC filter is applied to the analog data 
with a 3 dB point at 100 Hz. 

A DAQ program with real-time analysis running on 
seven PCs records candidate neutrino events using a 
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FIG. 1: Layout of the SAUND II array. Hydrophone locations 
are marked by bullets. Color shows the depth of the ocean 
floor. 



matched filter with the bipolar response function shown 
in Fig. [2] Because the phase response of the electronic 
system is unknown, SAUND II only considers the abso- 
lute value of the matched filter. Once per minute, the 
trigger threshold on each hydrophone is adjusted inde- 
pendently to a value that would have acquired 20 trig- 
gers/min in the previous minute. This differs from the 
SAUND I algorithm which only allowed thresholds to 
move up or down by a fixed step size. Each hydrophone 
triggers independently, and a 1 ms time series of only 
the triggered hydrophone is recorded. In addition, every 
minute the threshold on each hydrophone is recorded. 
The noise conditions in the form of a power spectral den- 
sity (PSD) and root mean square (RMS) are recorded ev- 
ery 5 s and 0.1 s respectively for each hydrophone. The 
SAUND II DAQ system consists of seven National In- 
struments cards (NI-PCI-MIO-64)and eight PCs of which 
seven are dedicated to data acquisition and one for con- 
trolling the entire system. It also includes a front-end 
module with hardware high-pass filters and an IRIG tim- 
ing signal distributor that connects to each of the cards. 
The NI ADC cards are controlled by COMEDI [1], and 
the DAQ software uses KiNOKO 0. 

By agreement with the US Navy, the SAUND II data 
acquisition system records data only when the array is 
not otherwise in use. From July 2006 to September 2007, 
the system ran under stable conditions for a total inte- 
grated time of more than 130 days as shown in Fig. [S] 

Thresholds used in SAUND II are compared to those 
used in SAUND I in Fig. 01 Since the same region of sea 
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FIG. 2: Bipolar response function ^ used for online trigger- 
ing. Markers show the digitized points used to match filter 
against the 156 kHz data sampling. 
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FIG. 3: Accumulated livetime of the SAUND H experiment. 



is measured over several seasons, ambient noise in both 
experiments is expected to produce similar distributions. 
However, SAUND II measures lower levels of adaptive 
thresholds overall. This is due to a better electrical noise 
performance of the new array, owing to the absence of 
long coaxial cables and digital transmission of signals on 
fibers. 



III. DATA ANALYSIS 

Tabic U shows the total number of triggers accumu- 
lated and the size of the remaining data set after selec- 
tion cuts are applied. Cuts 1 through 4 are applied on 
single-hydrophone triggers. After the source position in 
the ocean is triangulated using time difference of arrival, 
characteristics of these multi-hydrophone acoustic events 
are used on further cuts to test their consistency with 
expected properties of neutrino showers. In order to es- 
tablish data analysis cuts that are minimally biased, an 
8% subset of data was initially used to establish the data 
reduction cuts. Because the noise environment can vary 
greatly depending on time of day and season, the subset 
was selected carefully to have a random distribution in 
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FIG. 4: Comparison of adaptive threshold levels used by 
SAUND I [H and SAUND II. Threshold values have been 
normalized to account for the different sampling frequency 
and digitization of the response function (Fig. [2} in the two 
experiments. The quieter conditions in SAUND II are due to 
the absence of long coaxial cables, and an improved adaptive 
threshold algorithm. 



time. Cuts 1 through 3 result from analysing the 8% data 
alone. Once the full data was processed, it was found 
that some triggers are artificially generated by cross-talk 
between read-out channels. This effect was not discov- 
ered in the reduced data set as it affects a very small 
fraction of triggers. Therefore, parameters for cuts 4 to 
7 have been modified on the full data set to eliminate 
triggers caused by cross-talk while maintaining efficiency 
for neutrino showers. The method used to validate this 
procedure is described in section ITVl 

In TablelH the number of Online Triggers is the total 
number of triggers recorded by the online DAQ system 
in all hydrophones. Minutes of individual hydrophones 
with high threshold (>0.04), high trigger rate (>500 trig- 
gers/min compared to the target 20 triggers/min set by 
the DAQ), high trigger rate on a nearest neighbor, or 
man-made noise in the water are excluded to select Qual- 
ity Triggers. Man-made noise can contaminate large 
parts of the array, and can be identified by narrow-band 
frequencies continuously dominating a wide area of hy- 
drophones. This condition occurred a total integrated 
time of 9 min. Furthermore, data periods where the in- 
put cables appear to be disconnected are also removed. 
Such periods have unphysical, low noise levels, and are 
attributed to Navy personnel disconnecting individual 
cables for diagnosing problems in their system. One hy- 
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Cut 


Triggers 


Events 


1. Online Triggers 


327.9M 


— 


2. Quality Triggers 


146.7M 


— 


3. Waveform Selection 


2,814,545 


— 


4. Single Phone Rate 


2,562,047 


— 


5. Triangulation 


6,605 


4,995 


6. Isolated Event 


1,227 


320 


7. Radiation Pattern 


8 


2 



TABLE I: Data reduction process. The remaining number of 
single-hydrophone triggers after each cut is shown in the sec- 
ond column. The third column shows the number of acoustic 
events which are formed by four triggers that triangulate to 
a point in the ocean. 



drophone had to be excluded entirely due to the majority 
of data belonging to this condition. For the remaining 48 
hydrophones, this condition excludes '^6% of integrated 
livetime. 

Waveform Selection is used to further reduce back- 
ground triggers. For each trigger, the online data acquisi- 
tion system records a 1 ms waveform centered around the 
time of trigger. The absolute value of the matched-filter 
time series, M(t), is calculated using the response func- 
tion shown in Fig. [2] for each time series. The Waveform 
Selection cut is based on two parameters from this time 
series: the peak area ratio (PAR) and the Gaussian width 
(GW). The PAR is obtained by dividing the value of M 
obtained at the time of trigger by the integral of M in the 
1 ms window; PAR:^ M (to) / J*°^°;^ ^^ M{t)dt where to 
is the time of trigger. This quantity measures the signif- 
icance of the trigger compared to its surrounding times. 
The GW is obtained by smoothing the time series of M{t) 
to obtain the envelope shape, fitting this to a Gaussian, 
and extracting its width. For neutrino-induced signals, 
the width should be ^40 fj,s. Many marine animals vocal- 
ize around frequencies of interest for this work. However, 
the signals originating from such sources tend to have a 
larger GW value. Fig. \5\ compares these parameters for 
two triggers from the data set. The distribution of these 
parameters and the cut region where neutrino showers 
are expected are shown in Fig. |6l Triggers consisting 
of at least one data sample that saturate the ADC at 
±3.15 Pa bypass this cut and are kept regardless of the 
PAR and GW parameters. 

The Single Phone Rate cut removes surviving trig- 
gers that occur on a single hydrophone clustered in time. 
Triggers accompanied by 10 or more triggers in any 5 s 
window are removed. In addition, triggers that occur 
less that 10 ms from each other are also removed. These 
rates are significantly higher than the 20 triggers/min 
target trigger rate set by the DAQ. Diffused, steady state 
UHE neutrinos are not expected to arrive and interact in 
bursts, so the removal of clustered triggers on single hy- 
drophones only affects the livetime. Marine animals and 
artificial sources are usually responsible for such trigger 
patterns. 
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FIG. 5: Examples of Waveform Selection. The top panels 
show the time series, P(t), around two different triggers from 
the SAUND II data set. The panels at the bottom repre- 
sent their matched filter time series M(t) calculated with the 
~76 /iS bipolar response function, shown in Fig. [21 in 6.41 /is 
time steps. The GW parameter is obtained by smoothing 
M(t) and fitting to a Gaussian. The trigger on the left yields 
PAR=0.13 and GW=36.6 /xs. The trigger on the right yields 
PAR=0.02 and GW=97.9 /is. Even though the trigger on the 
right has a higher matched-filter value at the time of trigger, 
the PAR and GW parameters correctly discriminate between 
a trigger consistent with a shower (left) and a background 
event (right). 
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FIG. 6: Distribution of the parameters and accepted region 
of the Waveform Selection (top panel). The peak corresponds 
to the Knudsen noise floor [^. The 8% subset of data used 
to establish this cut is plotted. The accepted region was op- 
timized by generating over 300,000 Monte Carlo signals from 
]^q22_]^q25 gY neutrino showers observed at angles -2°-2° at 
various distances and adding recorded noise (bottom panel). 



With the remaining triggers Triangulation is per- 
formed. The 6,605 triggers produce 4,995 source loca- 
tions in the ocean by requiring four hydrophones to have 
the appropriate difference in arrival times. The time dif- 
ference of arrival method used is described in [l[ and in- 
cludes the effects of the depth-dependent speed of sound 
in the ocean. Only hydrophones that are vifithin 8 km 
of the hydrophone of earliest arrival time are considered 
for each event. For each event, a local rectangular co- 
ordinate system tangent to the Earth surface, centered 
around the earliest arrival-time hydrophone is used. High 
rates of trigger on a single hydrophone in the coincidence 
time window ('-~^8 s) produce combinatorics such that a 
single trigger can be involved in multiple events. How- 
ever, by further excluding events that trigger the same 
hydrophone within a minute of each other, only 320 time 
Isolated Events survive. 

Since acoustic radiation from neutrino showers is not 
expected to be emitted spherically, but rather in a disk- 
like shape orthogonal to the direction of the shower [2, la] , 
only showers of certain energy and orientation combina- 



tions can produce pulses capable of triggering all four hy- 
drophones from the triangulated event location. This ge- 
ometric constraint utilises the thresholds set at the time 
of the event, and only selects showers that are capable 
of triggering all four hydrophones with its acoustic radi- 
ation. However, not only does the radiation need to ex- 
ceed the thresholds, but estimated pressure amplitudes 
must also be observed in the recorded waveforms. Only 
two events have observed pressures matching or exceed- 
ing those expected by geometrically possible showers, and 
pass the Radiation Pattern cut. A metric, Rmissmg, 
quantifies the ratio of pressure missing in the recorded 
peak pressure, Pdet, from the estimated peak amplitudes. 

Pest- 



R 



missing 



T^UiPestf 



(1) 



where 



AP 



Pr 



det 



for Pest > Pdet 
for Pest < Pdet 



(2) 



Since AP is defined as the pressure missing, it is set 
to zero for waveforms meeting or exceeding estimated 
pressure amphtudes. R 



■raissing 



is calculated for a dense 
grid of shower energy-orientation combinations for each 
event. Fig. [7] shows the minimum values obtained by 
considering all possible Esh^ 0, and (j) configurations for 
each of the 320 events. 244 events did not have pos- 
sible showers with any direction and energy less than 
5x10^^ eV that fit the geometric constraint, and are as- 
signed Rmissing = 1-0. EvCUtS with Rraissing < 0.02 paSS 
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FIG. 7: Error metric Rmissing- The overflow bin at 
Rmissing = 1.0 coutains the 244 events with no possible 
shower geometries. The Monte Carlo events in the low- 
est bin extend beyond the maximum entries plotted to 131 
events. The inset shows details of the distribution close to 
zero. Again, the number of Monte Carlo entries in the lowest 
bin extends to 128 entries. 

the Radiation Pattern cut. The two events that pass re- 
quire the shower energy to be 10^** eV<£'s/i<5xlO^'' eV 
and 10^^ eV<£'s;i<5xlO^^ eV in order to be consistent 
with measured thresholds and peak pressures. 



IV. EFFICIENCY AND SENSITIVITY 
ESTIMATES 

The SAUND II detector is large enough to observe very 
different noise environments in different parts of the ar- 
ray. Therefore, the Quality Triggers cut described above 
removes livetime from single hydrophones minute-by- 
minute. Hence the configuration of enabled hydrophones 
is a complex pattern that varies from one minute to the 
next. To properly account for this, and to estimate the 
efficiency of the analysis cuts described above, neutrino 
events of different energies are generated by Monte Carlo 
simulation and injected randomly into the SAUND II 
livetime. 

The Monte Carlo study uses '^66M simulated neutrino 
events generated with the following distribution in zenith 
angle 9 and depth z. 



where I is the energy-dependent interaction length given 
by / = 1/apNA- o is assumed to scale with energy as 
CT ^ (Ey{GeV^^-'^^'^ [8], p is the density of sea water, 
and N A is Avogadro's number. The events are generated 
uniformly in the azimuthal angle as well as in latitude 
and longitude, in a 35 kmx57 km rectangular area com- 
pletely encompassing the SAUND II array. The closest 
hydrophone-boundary distance is >3 km. The fraction 
of neutrino energy deposited into the hadron shower, y, 
is generated with the distribution 



ji{y) = da /dy total = da/dycc + da/dy„ 



(4) 



/i 



Ae-'/^'""'^sin0 



(3) 



where da/dycc = (3 H- 2(1 — y)^)j/~°'®^ and da/dync = 
(1 -I- (1 — y)^)?/^"'®^ [9]. The hadron shower energy is 
obtained by Egh — y 'x E,^. We assume the same dis- 
tributions for all flavors of neutrinos. Ray tracing is 
performed between each hydrophone-event pair includ- 
ing the effects of sound speed variation in the ocean to 
obtain arrival times. The topography of the ocean floor 
is also considered. This produces some volume of ocean 
from where no rays can reach the hydrophones without 
intersecting with the ocean floor. An acoustic pulse is 
generated for each hydrophone-event pair using shower 
parameters [lO| and attenuation. For each event, a time 
of occurrence is assigned with a distribution that is uni- 
formly random in the SAUND II livetime. This allows 
the assignment of realistic noise conditions from the data 
to each hydrophone in each event. 

Trigger candidates are then passed through the cuts 
in Table ID First the matched filter values of the Monte 
Carlo candidate triggers are compared to the threshold 
set at the occurrence time and hydrophone. The Qual- 
ity Triggers cut keeps only the triggers that were as- 
signed a timestamp with good run conditions, while the 
Waveform Selection cut discriminates on simulated pulse 
shapes with added noise. 

The Single Phone Rate cut which is based on the rate of 
triggers on a single hydrophone can be affected by cross- 
talk. When a triggered waveform on one hydrophone 
produces cross-talk on other hydrophones that are also 
triggered, additional artificial triggers are created. The 
false increase in the trigger rate during this short period 
of time can affect the Single Phone Rate cut. To account 
for these situations, cross-talk triggers are generated by 
using the arrival times of surviving triggers and creating 
additional triggers on hydrophones within 8 km (elec- 
tronic cross-talk is instantaneous compared to both the 
speed of sound in water and the sampling frequency of 
the ADC cards). The Single Phone Rate cut and Tri- 
angulation are then performed including these false trig- 
gers. None of the artificial cross-talk triggers were found 
to generate events. The parameters used for the Single 
Phone Rate cut were chosen to minimize the rejection of 
true events in the presence of cross-talk. 

The Isolated Event cut creates deadtime that is much 
less than 1% of the integrated livetime. Finally Radiation 
Pattern fitting is performed to confirm that the metric 



R-n 



< 0.02, as shown in Fig. [71 



The efficiency of the SAUND II detector, 77 = 
ridetectediE) /ugeneratediE), is obtamcd at CEch dccadc of 
neutrino energy and shower energy. The effective vol- 
ume, Veff = Vgen X Vj wherc Vgen is the volume in which 
Monte Carlo events are generated, is shown in Fi g. |8l 
Using the model- independent method described in [ll|, 
the sensitivity for UHE neutrino flux at different neutrino 
and hadron shower energies is plotted in Fig. [S] using a 
90% confidence interval and assuming events remaining 
after cuts. 
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FIG. 8: The effective volume of the SAUND II experiment 
for different neutrino and shower energies. 



V. DISCUSSION AND RESULTS 

In order to formalize the SAUND II results, the two 
events surviving all cuts must be accounted for. The four 
triggers associated with each event are shown in Fig. [TU] 
and Fig. 1111 Further analysis capable of confirming the 
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FIG. 10: Waveforms of the four triggers associated with an 
event occurring 29 September 2006, 22:36:01.62 UTC com- 
patible with shower energy 5x10 eV>-Esh>10 eV, zenith 
angle 12.0°<6l<16.0°, and azimuth 170.7° <(/)<189.1° north. 
Three of the triggers saturate the ADC at ±3.15 Pa. Be- 
cause the phase response of the electronic system is unknown, 
SAUND II only considers the absolute value of the matched 
filter, and the phase of the bipolar pulse is not used in this 
analysis. 
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FIG. 9: Sensitivity to UHE neutrino flux at 90% confidence 
level obtained by Monte Carlo study of the SAUND II ex- 
periment assuming no events are detected. The sensitivity is 
plotted for different neutrino and hadronic shower energies. 
Although on average neutrinos create hadronic showers that 
have ~20% of its energy, the two plots are not simply a 20% 
shift from each other. This is due to the efliciency not scaling 
linearly with energy. 



FIG. 11: Waveforms of the four triggers associated with an 
event occurring 3 May 2007, 15:01:04.17 UTC compatible 
with shower energy 10^^ eV<ii'sh<5xlO^^ eV, zenith angle 
6.3°<6I<6.9°, and azimuth 233.8°<<7!)<235.5° north. 

remaining two events as background or signal will require 
better understanding of the detector. Because the phase 
response of the system is unknown, the phase informa- 
tion of the triggered pulses is not utilised to discriminate 



events. Thus the event shown in Fig. [TOl contains bipolar 
pulses with flipped phases, which are probably inconsis- 
tent as waveforms due to the same acoustic source. Apart 
from the phase, three of the four waveforms in Fig. 1101 
have similar high-amplitude bipolar pulses followed by a 
few more cycles of oscillation, while the fourth waveform 
on the top right has a much lower signal-to-noise. Since 
the probability of triangulating random arrival times to 
a single source location is extremely low, and since the 
four triggers are of order ~100 ms to seconds apart, the 
origin of this event is assumed to be acoustic despite its 
peculiar features. Furthermore, the three high-amplitude 
waveforms are similar to a type of background repeatedly 
seen in the SAUND II data set. All but the waveforms 
in Fig. [10] are eliminated due to their high repetition rate 
during their occasional bursts. The triggers in Fig. I10[ 
however, occur isolated in time and cannot be eliminated 
in such way. 

As discussed in [12], the phase response of a hy- 
drophone read-out system can distort bipolar pulses into 
more cycles, creating multi-polar oscillation signals. For 
this reason, the Waveform Selection cut does not define 
or utilize the number of cycles of oscillation. Thus some 
of the waveforms shown in Fig. llll are not clearly bipolar. 
In addition, because many triggers occur with amplitudes 
that are not significantly higher than the noise level, it 
is difficult to define how many peaks constitute a sig- 
nal. This effect also compounds the difficulty in phase 
matching bipolar pulses. 

Given the a posteriori observation made on two events, 
and the inability to perform a meaningful statistical anal- 
ysis with such a small data set, an upper limit on the flux 
is set. Since the knowledge of impulsive backgrounds is 
incomplete, an upper limit with statistical significance 
on neutrino flux alone cannot be derived. Nevertheless, 
a conservative limit can be obtained by assuming that the 
two events are due to signal, and no other background is 
present. Following the treatment in [lli] . 



\{e)^{e)de 



(5) 



where s is the expected number of events, A is the sen- 
sitivity of the experiment at different energies, and $ is 
the UHE neutrino flux. For a 90% confidence level. 



/ \{e)^{e)de < Sup 



(6) 



where Sup is the 90% Poisson confidence interval up- 
per limit of n event detection. At each neutrino energy 
bin, Sup{Esh) is calculated for different shower energy 
ranges Esh- For shower energy ranges lO^^'-lO^^ eV and 
10^^-10^'^ eV, n = 1 is used, making the upper bound 
Sup = 3.9. For other shower energy ranges, n = is 
used, making the upper bound Sup = 2.3. At each neu- 
trino energy, a weighted average (s„p) is calculated using 



{Sup{Eu)) = 2_^ Sup{Esh) 



NjEsh) 
N{E,) 



(7) 



where N{Esh) is the number of neutrinos that will pro- 
duce a shower in the energy range Eg^ when N{E^) neu- 
trinos are generated at neutrino-energy range E^. The 
weighted average gives {sup{E^)) at each neutrino energy, 
and is used to set a flux limit 



HE, 
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(8) 



The limit obtained by this method along with other ex- 
perimental limits are plotted in Fig. [121 
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FIG. 12: Neutrino flux upper limit from the SAUND II 
experiment. Various other limits are plotted: SAUND 
I H, GLUE [Hi, forte [H, ANITA [ij, ATCA- 
LUNASKA [III, and NuMoon-WSRT ^,^. 



VI. CONCLUSION 

The main factors contributing to the difference in sen- 
sitivities of SAUND I and SAUND II are volume, thresh- 
old, and hydrophone spacing. While SAUND II is advan- 
tageous in the first two factors, its sensitivity is limited 
by the larger hydrophone spacing. In both arrays, hy- 
drophones are mounted on the ocean floor, resulting in 
a mostly flat, two dimensional geometry which severely 
limits the sensitivity to the emission properties of show- 
ers, as already discussed in [l|. 

Unambiguously identifying rare UHE neutrinos from 
either acoustic or radio signals in huge naturally occur- 
ring bodies subject to varying and poorly understood 
backgrounds is very difficult. Combining multiple tech- 
niques in a hybrid array featuring acoustic, radio, and/or 
optical detectors makes it possible to identify neutrinos 
by detecting two signals that have very different produc- 
tion and propagation mechanisms as well as backgrounds. 



One possible array has been simulated in [1^ . In the ab- 
sence of such arrays, it is important to understand and 
identify backgrounds. The extensive study of the am- 
bient noise condition using SAUND II data ^ revealed 
that deep oceans are quieter than usually assumed at the 
high acoustic frequencies of interest. A model of such 
noise and a parametrization useful for ultra-high energy 
neutrino detection was developed. In executing a full 
analysis capable of detecting very rare events consistent 
with shower-induced signals, it was found that transient 
noise originating from discrete acoustic sources in the 
water causes the main limitation. A powerful cut, the 
Radiation Pattern cut, was developed to distinguish the 
geometry of emission to cope with such transient events. 
This cut will be significantly more effective when used 
with arrays of hydrophones arranged in all three dimen- 
sions. From such arrays, more information regarding the 



geometry of the emission lobe can be extracted. 
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